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ABSTRACT 


An analytical method is presented whereby an electrical resistance, flowing gas 
heater can be designed using wire mesh as the heating elements. By specifying the gas, 
the flow rate, the system pressure, the desired inlet and outlet gas temperatures, and 
the power supply characteristics, a compatible wire mesh heater element design can be 
determined. Hie calculations involve solving simultaneous equations for the heat gener- 
ation rate, the convective heat transfer from the mesh to the gas, and the change in en- 
thalpy of the gas at every axial position along the heater length. The method also enables 
the determination of the operating conditions of an existing mesh heater when used at some 
off-design condition or with another gas. 
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ANALYTICAL INVESTIGATION OF ELECTRICALLY HEATED WIRE MESH 

FOR FLOWING GAS HEATERS 
by William F. Mattson and William L. Maag 
Lewis Research Center 

SUMMARY 

An analytical method is presented whereby an electrical resistance flowing gas heater 
can be designed using wire mesh as the heating elements. By specifying the gas, the flow 
rate, the system pressure, the desired inlet and outlet gas temperatures, and the power 
supply characteristics, a compatible wire mesh heating element design can be determined. 
The calculations involve solving simultaneous equations for the heat generation rate, the 
convective heat transfer from the mesh to the gas, and the change in enthalpy of the gas at 
every axial position along the heater length. By varying the wire diameter, pitch, and 
size of the wire mesh, the designer is able to match the heater conditions to practically 
any type of power supply. The analytical method permits optimizing the heater design for 
compactness for specified operating conditions and limiting conditions such as maximum 
surface temperature or heat flux. It also enables one to determine the operating condi- 
tions of an existing mesh heater when used at an off- design condition or with another gas. 


INTRODUCTION 

High temperature, out-of-pile tests of nuclear fuel elements have been conducted at 
the Lewis Research Center using hot flowing gas to simulate reactor operating conditions. 
Electrical resistance heaters were developed to heat the flowing gas. One type of heater 
consisted of electrically heated wire mesh through which the gas flows. The analytical 
investigation leading to the design of the wire mesh heating elements is the subject of this 
report. 

The specific requirement was to design a heater utilizing an existing power supply to 
heat hydrogen, helium, and nitrogen gases at specified flow rates to temperatures up to 
2780° K. The search for a suitable heating element led to a commercially available tung- 
sten wire mesh made of helical coils. The mesh provided the capability of obtaining 



electrical resistance compatible with the current and voltage ratings of available power 
supplies and also possessed a large surface area for heat transfer. Tungsten with its 
high melting point was chosen because of the high gas temperature requirement. A liter- 
ature survey revealed only a limited amount of low temperature experimental heat trans- 
fer data for wire mesh. References 1 and 2 report the results of a Lewis experimental 
program to establish high temperature forced convection heat transfer correlations for 
electrically heated wire mesh. 

After a heat transfer correlation was established, the task of using the mesh to de- 
sign a heater capable of delivering gas at the desired flow rate and temperature remained. 
This involved solving simultaneous equations for the heat generation rate, the convective 
heat transfer from the mesh to the gas, and the change in enthalpy of the gas at every 
axial position along the heater length. This report investigates the variables in mesh 
geometry that can be utilized to achieve compact heater designs for limiting conditions 
such as maximum mesh surface temperature or heat flux. It defines reasonable limits 
within which a mesh heating element can be expected to operate satisfactorily. Two power 
supplies representing high and low electrical resistances are considered to illustrate how 
a mesh can be adapted to each. Experimental data from an operating mesh heater are 
presented to corroborate the analytical results. 

A computer program by Geraldine E. Amling is included as appendix B. 


METHOD OF ANALYSIS 

In designing an electrical resistance gas heater certain operating conditions must be 
specified. These include the type of gas, the gas flow rate and pressure, the initial and 
final gas temperature, and the current and voltage rating of the power supply. The next 
step then is to design a heating element that can generate and transfer heat to the flowing 
gas at a reasonable heat flux and at a surface temperature that would ensure a reasonable 
heater life. For this investigation, two different power supplies are considered. One 
represents a high resistance system in which the current can be maintained between 150 
and 1500 amperes at voltages up to 3750 volts. The second represents a low resistance 
system whereby the voltage can be controlled between 10 and 50 volts at currents up to 
60 000 amperes. The working gas is hydrogen at a flow rate of 0. 0453 kilogram per 
second with an inlet temperature of 278° K and a desired outlet temperature of 2780° K. 


Mesh Heating Element 

Tungsten wire mesh that can be electrically heated to surface temperatures greater 
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Figure l. - Mesh heating element. 


than 2780° K becomes a unique heater element when a nonelectrically conductive fluid 

flows through it. The mesh provides a large heat transfer surface area per unit volume 

of heater element that results in high power densities and small compact heaters. The 

7. 6- by 2. 5- by 0. 36-centimeter mesh shown in figure 1 has a surface area to volume 
2 3 

ratio of 19 (cm /cm ). The surface area of this mesh is equivalent to that of a 2. 54- 
centimeter-diameter tube, 16. 3 centimeters long. The chief disadvantage of a mesh ele- 
ment, as with any porous media, is the pressure drop restriction it places on fluid flowing 
through it. If relatively high pressure drops can be tolerated, then wire mesh provides 
many of the features necessary in a good heater element. 

The wire mesh heater element considered herein is represented by figure 1. The 
element consists of helical coils of tungsten wire interwound into a mesh. Both ends of 
the coils are sandwiched between two tungsten plates each 0. 153 centimeter thick. The 
plates and the coil ends are heliarc welded together at both ends of the mesh. The plates 
provide electrical bus connections for directing current through the coils. The total cur- 
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rent to the end plates will distribute itself into each wire coil according to the resistance 
of that coil. If the coolant flow distribution is uniform throughout the mesh, then the heat 
flux will be uniform and each coil will carry an equal amount of current. This is a neces- 
sary assumption for the analysis made herein although experimental results indicate that 
such is not always true. Reference 2 states that interwound mesh heating elements are 
susceptible to hot spotting and subsequent burnout at conditions of high heat flux with a 
nondissociating gas such as helium. This is believed to be caused by slight nonuniform- 
ities in the interwound geometry resulting in flow maldistribution within the mesh. If the 
mesh heat flux is maintained below some predetermined critical value, then the previous 
assumption is correct. 

In order to calculate the heat transfer from a mesh it is necessary to know such 
geometrical parameters as surface area, current cross-sectional area, current path 
length, porosity, and equivalent diameter. These will be defined in terms of wire diam- 
eter, inside coil diameter (commonly referred to as mandrel diameter), coil pitch, num- 
ber of parallel coils per width of mesh, and mesh length. The length S of wire in a 
single helical coil is given in reference 3 as 


S = - y ct 2 (D + d) 2 + p 2 (1) 

where b is the mesh length between end plates, p is the coil pitch, D is the inside coil 
diameter, and d is the wire diameter. (All symbols are defined in appendix A. ) The 
total heat transfer surface area for N parallel coils is 

A g = ?rdSN (2) 

This was the area used in determining the heat transfer coefficients of references 1 
and 2. The total current cross-sectional area for N parallel coils is 

A = 1^? N (3) 

c 4 


The average mesh porosity is defined as the ratio of void volume to total volume or 


e 


bw(D + 2d) - — SN 
4 

bw(D + 2d) 


( 4 ) 
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where w is the mesh width and (D + 2d), a coil diameter, respresents the thickness in 
the flow direction. The equivalent diameter for a porous structure is usually defined as 


£ _ 4(void volume) _ 4bw(D + 2d)e ^ 

e surface area A 

S 

The mass velocity of fluid through the element is based on an average flow area within the 
mesh as defined by the product of frontal area and porosity; that is, 


G = 


W 

bwe 


( 6 ) 


For this analysis the mesh completely fills the flow passage so that there are no fluid by- 
pass effects as discussed in reference 1. 

For this study, certain ground rules were established so that the mesh designs used 
were as compact as possible within reasonable limits for fabrication and assembly. For 
this reason the interwound helically coiled mesh of reference 1 was specified. The inside 
diameter of the coils was defined as 2d + 0. 005 centimeter. The 0. 005 centimeter allowed 
clearance for interwinding the coils. The coil pitch was defined as 3d which represents 
the minimum allowable pitch for interwinding. The number of coils N can be expressed 
in terms of mesh width and the wire and mandrel diameters according to the following 
empirically derived equation: 

N = 3^-- (2d + °>1 (7) 

L 2d + D J 

This relation states that for a width w, less a half coil width on each side, the number 
of interwound parallel coils is three times the number that would be available if the coils 
were just placed side by side. Equation (7) is an approximation which is only to be used 
when coil pitch and mandrel diameter are defined as done previously. It is necessary 
that N always be evaluated as a whole number for any actual design. 

With the foregoing set of rules different meshes having the same frontal dimensions 
were developed merely by changing the wire diameter, the other parameters all being 
defined in terms of this variable. Wire diameters from 0. 036 to 0. 127 centimeter were 
used for this study. A mesh made from wires smaller than 0. 036 centimeter in diam- 
eter was considered from previous operating experience to be too flimsy, while a mesh 
made from wires having a diameter greater than 0. 127 centimeter was considered too 
rigid for tungsten fabrication. The geometrical properties of representative meshes are 
given in table I for a 6. 5-square-centimeter mesh (2. 54 cm wide by 2. 54 cm long). Prop- 
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TABLE I. - MESH GEOMETRY PARAMETERS 


Wire di- 

Mandrel 

Number of 

Coil 

Mesh 

Ratio of total 

Ratio of total 

Equiva- 

ameter, 

diameter, 

parallel coils, 

pitch, 

porosity, 

surface to 

wire length to 

lent di- 

d, 

D, 

N, 

P, 

€ 

frontal area, 

current cross- 

ameter, 

cm 

cm 

coils/2. 54 cm 

cm/turn 


A s /Af, 

sectional area, 

D e> 






cm 2 /cm 2 

S/A c’ 2 

cm/cm 

cm 

0.036 

0.076 

49 

0. 107 

0.56 

7. 37 

180 

0.045 

.051 

. 107 

34 

. 152 

. 56 

7. 17 

127 

.066 

.076 

. 157 

22 

.228 

. 58 

6.84 

87 

. 105 

. 102 

. 208 

16 

. 304 

.60 

6. 52 

68 

. 151 

. 117 

. 239 

13 

. 350 

.61 

6. 34 

60 

. 181 

. 127 

. 259 

12 

.381 

.62 

6. 21 

56 

. 204 


erties such as surface area, equivalent diameter, and porosity relate to the heat transfer 
characteristics of the element while S/A c is a measure of the electrical resistance. 

Mathematical Model 

A physical model of two typical mesh stages is shown in figure 2. The gas enters 



Figure 2. - Physical model. 
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upstream of stage 1 at temperature T^. On passing through the electrically heated mesh 
of stage 1, heat is transferred to the gas by forced convection and exits at Tg; that is, 

Q = hA g (T s - T b ) (8) 


where 


T b = 


T l + T 2 


( 9 ) 


The heat transfer coefficient h for interwound mesh was experimentally determined in 
reference 1 to be of the following form: 


h 



Substituting equations (9) and (10) into equation (8) gives 


Q = 



T 1 + T 2j 


( 10 ) 


( 11 ) 


Viscosity fi and thermal conductivity k values were obtained as functions of tempera- 
ture and pressure from reference 4. Specific heat at constant pressure was evaluated at 
the stage inlet gas temperature using the specific heat-temperature correlation of refer- 
ence 5. The electrical power input to the mesh heating element is given by the following 
equation in which either the voltage or current is specified depending upon the type of 
power supply available: 


Q =r =I 2 R 

R 


where 



( 12 ) 


(13) 
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Resistivity p for most metals can usually be expressed as some function of temperature 
(ref. 6). Substituting equation (13) into equation (12) gives 

V 2 

Q = — — (14) 

S 


The enthalpy rate increase of the gas is given by 

Q=WC p (T 2“ T l ) (15) 

The simultaneous solution of equations (11), (14), and (15) for a given set of conditions 
establishes the steady state heat transfer for the mesh- gas -power system. The three 
unknowns in these equations are the surface temperature T , the outlet gas temperature 
from stage 1 Tg, and the rate of heat transfer Q. The numerical method of reference 7, 
chosen to solve these equations, was adapted to a digital computer. Appendix B presents 
details of the method as well as the computer program used. This analysis assumes an 
adiabatic system. 

The computer program solves these equations for stage 1 and then proceeds to the 
next stage using the outlet gas temperature from stage 1 as the new inlet temperature. 

The calculations proceed stage by stage until the desired gas temperature is obtained or 
until certain predetermined limits such as surface temperature, heat flux, or number of 
stages are exceeded. When limits are reached, it is necessary to revise the mesh ele- 
ment design and repeat the calculations. The variables in mesh geometry that affect the 
calculations are discussed in the following section. 


RESULTS AND DISCUSSION 


The most compact gas heater results when the heat transfer surface operates at a 
constant, maximum temperature. For a tube heater the equation 



(16) 


gives the minimum length to diameter ratio for the case of constant wall temperature as 
derived by Reynold's analogy (ref. 8). For the design flow rates of this report, a tube 
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friction factor of 0. 005 is reasonable. This results in an L/D^. of about 200. To elec- 
trically heat a tube of this length and maintain a constant wall temperature requires that 
the current cross-sectional area vary throughout the tube length. This is very difficult 
and impractical to do. 

The advantages of a mesh heating element are compactness and design flexibility. 

The latter permits the designer to more or less specify the heat generation rate or the 
surface temperature at any axial position along the heater length. This may be accom- 
plished by dividing the heater length into mesh stages in which the heat generation, heat 
transfer, and enthalpy for each stage or group of stages is specified. Then, by varying 
the mesh surface area, equivalent diameter, and ratio of current path length to cross- 
sectional area, it is possible to approximate the desired surface temperature or heat flux 
conditions. Referring to table I reveals that the surface area of the meshes investigated 
varied about 20 percent while the equivalent diameter and S/A varied by factors of 4. 6 
and 3. 2, respectively. 

Two design cases were investigated. One required a high electrical resistance mesh 
heater element, while the second required a low resistance element. Both cases specified 
hydrogen gas at 0. 045-kilogram-per- second flow rate with an inlet temperature of 278° K 
and a desired outlet temperature of 2780° K. The system pressure was 100 atmospheres 
so that hydrogen dissociation was not a factor. The heater flow area was 25. 8 square 
centimeters, and the mesh stage completely filled the flow passage cross-sectional area. 
The surface temperature of the mesh could not exceed 3200° K, and no more than 20 mesh 
stages were used. There was no limitation put on heat flux, but references 1 and 2 show 
that their maximum experimental value was about 1. 24 kilowatts per square centimeter 
and the operating experience of reference 9 indicates that an average heat flux of 0.62 kilo- 
watt per square centimeter seems reasonable for an extended period of operation. The 
heat transfer correlation of reference 1 was used for the calculations. 


High Voltage Heater Design 

For a high voltage power supply the heater elements are connected electrically in 
series such that the same current flows in each mesh. The heat generated in each mesh 
stage is therefore proportional to the S/A c ratio and the resistivity of tungsten at the 
temperature of the element. Figure 3 represents this type of heater. For this particu- 
lar design each mesh stage is segmented into four identical panels so that a high electri- 
cal resistance can be obtained. The current flows through each of the four panels in 
series before it passes onto the next mesh stage. Electrical insulating material (boron 
nitride in this case) separates the panels so that each channel directs one-fourth of the 
gas flow axially through the heater. Since each mesh panel in a stage is identical and the 
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Figure 3. - Series heater. 




channel flow areas are equal, the heat flux may be assumed equal throughout the stage. 
The mesh dimensions used were 1. 27 centimeters wide by 5. 08 centimeters long for each 
of the four flow passages in one stage for a total mesh frontal area of 25. 8 square centi- 
meters per stage. 

Figure 4 shows the effect of wire diameter on heater design. The large diameter 
(0. 127 cm) wire mesh has such a low electrical resistance that the heat flux is very 
small. A heater using this size wire throughout would be prohibitively long to achieve 
the desired gas temperature. The small diameter (0. 036 cm) wire mesh has a much 
higher electrical resistance with about the same surface area. Consequently, the mesh 
operates at a higher surface temperature that increases exponentially with heater length. 
This is caused chiefly by the increase in tungsten resistivity with temperature. The re- 
sult is a short compact heater in which the major portion of the heat is generated and 
transferred in the last few stages. This is not a good design because it is susceptible to 
mesh burnout. The combination of high surface temperature, high heat flux and small, 
hence weak, wires would probably result in a burnout in the event of any small maldistri- 
bution in gas flow through the mesh. 

Figure 5 shows the effect of maximizing surface temperature. This was done by 
choosing mesh from table I that would give the highest electrical resistance for each 
stage without exceeding 3200° K. The result is a very short, compact heater in which 
the heat flux is very high. Since these heat fluxes have not been demonstrated experi- 
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Voltage, V Heat flux, kW/cm‘ Temperature, 




(a-1) Wire diameter, 0.036 centimeter. 





(b-1) Wire diameter, 0. 036 centimeter. 
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(a-2) Wire diameter, 0. 127 centimeter, 
(a) Temperature profile. 
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(b-2) Wire diameter, 0. 127 centimeter, 
(b) Heat flux profile. 
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Heating stage number 


(c-1) Wire diameter, 0. 036 centimeter. 


(c-2) Wire diameter, 0. 127 centimeter, 
(c) Voltage profile. 


Figure 4. - Effect of wire diameter for high voltage series connected heater. Operating conditions: hydrogen gas ; 
gas flow rate, 0.045 kilogram per second; current, 1000 amperes. 
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(a) Temperature profile. 




Figure 5. - Effect of maximizing surface 
temperature for high voltage series 
connected heater. Operating condi- 
tions: hydrogen gas ; gas flow rate, 
0.045 kilogram per second; current, 
1500 amperes. 
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mentally, this does not represent a practical design. Another factor that would tend to 
negate this design would be the high voltage drops across each mesh. For example, the 
maximum voltage difference between adjacent stages of figure 5 is 682 volts. Providing 
sufficient insulation to prevent arcing could be a limiting problem, particularly at the 
peak temperatures. 

From the foregoing discussion it may be concluded that a practical series heater 
design would consist of various meshes that would all operate at a fairly uniform heat 
flux. The heat flux would be high enough so that the heater dimensions will remain com- 
pact but not high enough to cause burnout due to flow maldistribution. The smaller wire 
diameter mesh would be used at the heater inlet, and as the operating temperature in- 
creases the wire size would increase so that the mesh operating at the highest tempera- 
ture is also the strongest. 


Low Voltage Heater Design 

For a constant low voltage, variable high current power supply the mesh heater ele- 
ments are connected electrically in parallel. Each mesh operates with the same voltage 
drop but the current through each varies inversely with the resistance. Consequently, 
the cooler elements carry more current and operate at a higher heat flux than the hotter 
ones, just the opposite of the series heater. 

Figure 6 represents this type of parallel heater. For this analysis the mesh panels 
are 5. 08 centimeters by 5. 08 centimeters. The power is introduced through water cooled 
copper tubes and buses. The heater elements are bolted onto the buses which form two 
sides of the square flow passage. Insulating material (boron nitride) forms the other two 
sides of the channel. The gas flows in series through each mesh, increasing in tempera- 
ture from one stage to the next. 

Figure 7 shows the effect of wire diameter on heater design. The low electrical re- 
sistance of the large diameter (0. 127 cm) wire mesh makes it more compatible with this 
type of power supply than the small wire diameter mesh. The result is a short heater 
operating at a high average heat flux with the 0. 12 7- centimeter wire as compared to a long 
heater operating at a relatively low average heat flux for the 0. 036-centimeter wire. In 
comparing these designs with the series heater it becomes evident that the temperature 
and heat flux gradients for the parallel heaters are not as severe as those of the series 
heaters. This is because the parallel electrical connection has an equalizing effect. 

Since the resistance of each mesh increases with temperature, the current will distribute 
itself inversely with temperature. Consequently the cooler mesh will carry more current 
and the hotter ones less for the same applied voltage, thereby flattening the mesh temper- 
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Current, I, A Heat flux, kW/cirr Temperature, 




(a-1) Wire diameter, 0. 036 centimeter. (a-2) Wire diameter, 0.127 centimeter, 

(al Temperature profile. 




(b-1) Wire diameter, 0. 036 centimeter. (b-2) Wire diameter, 0.127 centimeter, 

(b) Heat flux profile. 




Heating stage number 

(c-l> Wire diameter, 0. 036 centimeter. (c-2) Wire diameter, 0.127 centimeter, 

(c) Current profile. 

Figure 7. - Effect of wire diameter for high current parallel connected heater. Operating conditions: 
hydrogen gas ; gas flow rate, 0.045 kilogram per second; voltage, 25 volts. 
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Figure 8. - Effect of maximizing surface 
temperature for high current parallel 
connected heater. Operating condi- 
tions: hydrogen gas; gas flow rate, 
0.045 kilogram per second; voltage, 

50 volts. 
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ature profile and to some extent the flux profile. Additional equalization can be achieved 
by varying the mesh wire diameter. 

Figure 8 shows the effect of maximizing surface temperature by tailoring each mesh 
stage. As with the series heater, the result is a very short compact heater with the 
mesh operating at very high heat fluxes. In addition to being susceptible to burnout, this 
design would impose very large demands on the electrical busing. The necessity of 
delivering 40 000 amperes to these five mesh stages in a small volume of space would be 
a very difficult design problem. 

The same general conclusion applies to the parallel heater as for the series heater. 
The better design is achieved by equalizing the mesh heat flux at some reasonable value 
and minimizing the number of mesh operating near the maximum allowable surface tem- 
perature. Reference 8 represents this type of compromise design. 


Effect of Gas 


A heater designed for one gas will operate quite differently with another gas. Fig- 
ure 9 compares nitrogen operation with the parallel heater designs of figure 7. The com- 
parison is made with a nitrogen flow rate that corresponds to the same power and gas 
temperature rise as for the hydrogen design conditions; that is, 



= W 


'P, H, 


H 


2 C 


p, No 


(17) 


For the same applied voltage, the mesh operate at a higher surface temperature for ni- 
trogen than for hydrogen. This is chiefly because nitrogen has a lower thermal conduc- 
tivity which causes its heat transfer coefficient to be less than that for hydrogen as seen 
by equation (10). By comparing thermal conductivities it is possible to predict heater be- 
havior for gases other than the design gas. As in the previous cases, it would be neces- 
sary to operate the nitrogen heater at both a lower power level and flow rate to achieve 
2780° K gas. However, in order to accurately predict the flow rate and power level, it 
would be necessary to use the computer program described in appendix B. 
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Current, I, A Heat flux, kW/cm^ Temperature, 



(a) Temperature profile. 



(b) Heat flux profile. 




(a-2) Wire diameter, 
0. 127 centimeter. 



(b-2) Wire diameter, 
0. 127 centimeter. 



(c-2) Wire diameter, 
0. 127 centimeter. 


(c) Current profile. 

Figure 9. - Effect of nitrogen gas on operation for high current parallel connected heater. 
Operating conditions: nitrogen gas ; gas flow rate, 0.63 kilogram per second; voltage, 

25 volts. 
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Experimental Results 


An electrical resistance wire mesh gas heater has been designed and built with the 
aid of the analytical program described herein. A detailed description of this heater is 
presented in reference 9. The heater was originally designed for hydrogen but has been 
used both as a nitrogen and hydrogen heater. It consists of four banks of mesh stages 
similar to the three banks shown in figure 6. Each bank represents a different mesh 
geometry and each is powered by a separate low voltage power supply. The mesh heating 
elements are connected in parallel as in figure 6. Figure 10 compares experimental re- 
sults for nitrogen operation with the calculated conditions from the analytical program. 
The calculated average surface temperatures were 1 to 9 percent higher than the experi- 
mental average values determined from resistivity measurements but the calculated and 
measured outlet gas temperature differed by less than 5 percent. These differences are 
within experimental error and do not detract from the significance of the analytical re- 
sults. 



Heating stage number 

Figure 10. - Comparison of calculated to experimental results for high 
current parallel connected heater. Operating conditions: nitrogen 
gas; gas flow rate, 0.313 kilogram per second. 
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SUMMARY OF RESULTS 


An analytical method is presented whereby an electrical resistance gas heater can be 
designed using wire mesh as the heating elements. By specifying the type of gas, the flow 
rate, the system pressure, the desired inlet and outlet temperatures, and a power supply 
of sufficient capacity and known voltage and current rating, a compatible wire mesh design 
can be determined. The calculations involve solving simultaneous equations for the heat 
generation rate, the convective heat transfer from the mesh to the gas, and the change in 
enthalpy of the gas at every axial position along the heater length. 

Wire mesh heating elements provide the designer with compactness and the flexibility 
to match the heater conditions to practically any type of power supply. For low voltage 
power supplies, the mesh are connected electrically in parallel while for high voltage 
power supplies, a series electrical connection is desirable. In both cases the mesh geom- 
etry parameters change with heater length to comply with the varying heat transfer condi- 
tions throughout the system. High heat fluxes and surface temperatures should be avoided 
since they increase the chance of mesh burnout because of slight flow maldistributions 
that are inherent with any porous geometry. 

Given specific operating conditions and limiting conditions such as surface tempera- 
ture and heat flux, the analytical method permits one to optimize the heater design for 
compactness. It also permits one to determine the operating conditions for the heater 
when it is used for a different purpose such as with another gas or at some off-design con- 
dition with the same gas. Comparing the results from the analysis with actual experimen- 
tal results of an operating heater demonstrates the usefulness of this analytical method. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 23, 1968, 

120-07-04-54-22. 
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APPENDIX A 


SYMBOLS 


A c 

total current cross-sectional area, 

o 

N 

number of parallel coils 


cm 

o 

n 

exponent in eq. (10) 

A f 

total frontal area, cm 

P 

coil pitch, cm/turn 

A s 

o 

total surface area, cm 

Q 

rate of heat transfer to gas, J/sec 

b 

length of coil, cm 

R 

electrical resistance, ohm 

C 

coefficient in eq. (10) 

S 

total wire length of helical coil, cm 

C P 

specific heat of gas at constant 

T b 

average bulk temperature, 
°K 


pressure, J/(kg)(°K) 


D 

mandrel diameter, cm 

T s 

average surface temperature, °K 

D e 

equivalent diameter, cm 

T w 

wall temperature, °K 

D t 

tube diameter, cm 

T 1 

stage inlet gas temperature, °K 

d 

wire diameter, cm 

T 2 

stage outlet gas temperature, °K 

f 

friction factor 

V 

voltage, V 

G 

average mass velocity, 
2 

(kg/sec)/cm 

W 

gas flow rate, kg/sec 

h 

average heat transfer coefficient, 

w 

mesh width, cm 


J/(sec)(cm 2 )(°K) 

e 

mesh porosity 

I 

current, A 

P 

absolute viscosity of gas, 

k 

thermal conductivity of gas, 


kg/(sec)(cm) 


J/(cm)(sec)(°K) 

P 

resistivity of tungsten, ohm-cm 


L tube length, cm 
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APPENDIX B 


COMPUTER PROGRAM 

by Geraldine E. Amling 


The computer program presented herein was written for the purpose of designing 
wire mesh heating elements for an electrical resistance flowing gas heater. If the gas, 
the flow rate, the system pressure, the desired inlet and outlet gas temperatures, and 
the power supply characteristics are known, a compatible wire mesh heater element can 
be determined. Basically, this program solves the following three simultaneous equa- 
tions having three unknowns, surface temperatures TS, outlet gas temperature T2, and 
the rate of heat transfer to gas Q: (English units, as defined in the section Computer 
Code Symbols of this appendix, are used throughout. ) 


Q =• 


0. 948x10" 3 E 2 


(Bl) 


5. 2X10 


•6/ TS 


v l. 185 


aoooj 


XLA • XL 


Q = — (ASF) 
144 


Q = FR* CP(T2 - Tl) 


0 364 l 144FR-DE \ 0 - 542 XK 
\P • AR • XMU/ DE 


TS - 


Tl + T2' 


(B2) 

(B3) 


If current is to be specified instead of voltage, equation (Bl) would be changed to 


Q = I 2 (5. 2x10” 6 ) 185 (XLA)(XL) (B4) 

\ 1000 j 

The resistivity of tungsten, 5. 2xl0"^(TS/1000)^‘ is from reference 6. Since the 
Prandtl number in reference 1 is nearly a constant, it was not included in the correlation 
presented in equation (B3). The method used to solve these equations was the Newton- 
Raphson iteration process (ref. 7) with Gauss double precision routine to solve the ma- 
trix. For the specific case presented, gas properties are evaluated at the mesh surface 
temperature; however, evaluation could be at any specified temperature with only minor 
changes to the program. The solution requires an initial estimate of TS, Q, and T2 in 
the region of convergence of the functions used. Corrections are then repeatedly calcu- 
lated and added to the values until the desired accuracy is obtained. If AH, AL, and 
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AK are corrections to the initial estimate for Q, T2, and TS, respectively, then 


Q1 + AH = Q’l 

(B5) 

T2 + AL = T’2 

(B6) 

TS + AK = T'S 

(B7) 


The new, more exact values are then used to calculate new corrections. 

In the program, preliminary estimates of starting values are obtained by using a 
single iterative process to determine, within 100° R, the value of TS, and the corre- 
sponding values of T2 and Q. The approximation is further refined using the subroutine 
MGAUSD to obtain a TS value which is within 1 percent of the true value. 

Starting values used in this program are approximations, and although satisfactory 
results were obtained in all cases run, there is a possibility that the Newton-Raphson 
method will not converge for a given set of values. 

Hydrogen, nitrogen, and helium properties have been curve fitted in the range of 
500° to 6000° R. The coefficients of the curves are included in the main program as 
data,, where the units of these properties are as follows: 

Temperature, °R; thermal conductivity XK, Btu/(ft)(sec)(°R); viscosity XMU, 
lb/(ft)(sec); specific heat CP, Btu/(lb)(°R). 

If the surface temperature TS becomes negative, a print out of TS is given, com- 
putation is terminated, and program reads in a new case. 

Accuracy of the approximations using Newton-Raphson technique was established by 
convergence of Q values at the 1. 0 percent level. However, the percent error may be 
increased or decreased according to the user's needs. Percent error tests are made at 
statement number 14 and statement number 16. If specified accuracy is not achieved 
after 50 iterations, a print out of TS, T2, and Q is given, computation is terminated, 
and program reads in a new case. 


Input 

The input format sheet is shown in table H. The quantities filled in are for a sample 
problem. The output for this sample problem will be presented subsequently. 
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TABLE n. - INPUT FORMAT FOR FORTRAN COMPUTER PROGRAM 





























Instructions for preparing the input are given in table HI. 


TABLE HI. - INSTRUCTIONS FOR PREPARING INPUT 


Card type 

Number 

Columns 

Description of variables 

1 

1 

1 to 2 

CODE 

two- letter code for type gas 





(H2, HE, or N2) 



3 to 4 

NN 

number of stage groups 

2 

NN 

1 to 8 

D 

wire diameter, in. 

(constants 


9 to 16 

E 

voltage, V 

for each 


17 to 24 

ASF 

A s /A f , total surface area/total 

stage) 




frontal area 



25 to 31 

DE 

equivalent diameter, ft 



33 to 40 

P 

mesh porosity, percent 



41 to 48 

XLA 

S/A c , total wire length/total 





current cross-sectional area 





(for 1 sq in. of mesh) 



49 to 56 

XI T 

last mesh number of each stage 

3 

1 

1 to 8 

PCE 

percent error 

(data for 


9 to 16 

AR 

total frontal area, in. ^ 

one 


17 to 24 

XL 

mesh length/mesh width 

case) 


25 to 32 

FR 

gas flow rate, lb/sec 



33 to 40 

Til 

T 1? stage inlet gas temperature, 





°R 


Card type number 2 may be repeated for a group of stages (NN). The sample prob- 
lem has a total of 30 stages. The first group contains 6 stages (XIT(l) = 6. 0); the second 
group, 8 stages (XIT(2) = 14. 0); the third group, 8 stages (XIT(3) = 22. 0); and the fourth 
group contains 8 stages (XIT(NN) = 30. 0). 

This program is written in FORTRAN IV and is operational on the IBM 7094-2/7044 
direct-coupled system of the Lewis Research Center. A flow chart of the program is 
shown in figure 11. 

Machine running time is approximately 2 minutes per case. 




Figure 11. - Flow chart for main program HEWMAN. 
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Output 


Input and output are printed out for each iteration. Output consists of the following: 
T2 outlet gas temperature, °R 
TS surface temperature, °R 
Q rate of heat transfer to gas, Btu/sec 
RT total resistance 
IT mesh number 


A1 

A2 

A3 

AH 

AK 

AL 

AR 

ASF 

CODE 

CP 

CPA 

CPAH 

CPAHE 

CPAN 

CPBPH 

CPBPHE 

CPBPN 

CTS 

D 


Computer Code Symbols 

first equation (Bl) 
second equation (B2) 
third equation (B3) 
correction factor for Q 
correction factor for T 

s 

correction factor for T2 

9 

Af, total frontal area, in. 

2 2 

A g /A f , total surface area/total frontal area, in. /in. 
two letter code for element (H 2 , He, or N 2 ) 
specific heat of gas at constant pressure, Btu/(lb)(°R) 
current coefficients, C curve, at Til 

ir 

coefficients, C curve, hydrogen 

r 

coefficients, C p curve, helium 
coefficients, C curve, nitrogen 

r 

breaking points, C curve, hydrogen 

r 

breaking points, C curve, helium 

r 

breaking points, C curve, nitrogen 

r 

current value of T 
s 

d, wire diameter, in. 
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DCP 

DE 

DK 

DM 

DRES 

E 

FR 

IT 

KKC 

L 

NN 

NNTEST 

P 

PCE 

Q 

QQQ 

RES 

RHO 

RT 

TB 

TS 

TS1 

Til 

T2 

T21 

XI T 

XK 

XKA 

XKAH 


derivative of C at Til 
P 

D e , equivalent diameter, ft 
derivative of k at CTS 
derivative of / u at CTS 

reciprocal of electrical resistance at each stage, ohm 
voltage, V 

W, gas flow rate, lb/sec 
XIT, mesh number 
stage number 
denotes gas used 
number of stage groups 
last stage of each case 
e , mesh porosity, percent 
percent error 

rate of heat transfer to gas, Btu/sec 
initialization of Q 

l/RT, reciprocal of electrical resistance 

p, resistivity of tungsten, ohm-ft 

total electrical resistance, ohm 

T^ = (Tj + Tg)^, average bulk temperature, °R 

T_, average surface temperature, °R 
s 

initialization of T 
s 

initial Tj, stage inlet gas temperature, °R 
stage outlet gas temperature, °R 
initialization of Tg 
mesh number 

k , thermal conductivity of gas, Btu/(ft)(sec)(°R) 
s 

current coefficients, k curve, at CTS 
coefficients, k curve, hydrogen 
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XKAHE 

XKAN 

XKBPH 

XKBPHE 

XKBPN 

XL 

XLA 

XM 

XMA 

XMAH 

XMAHE 

XMAN 

XMBPH 

XMBPHE 

XMBPN 

XMU 


coefficients, k curve, helium 
coefficients, k curve, nitrogen 
break points, k curve, hydrogen 
break points, k curve, helium 
break points, k curve, nitrogen 
mesh length/mesh width 

S/A total wire length/total cross-sectional area (for 1 sq in. of mesh) 
p at CTS 

current coefficients, p curve, at CTS 
coefficients, p curve, hydrogen 
coefficients, p curve, helium 
coefficients, p curve, nitrogen 
break points, p curve, hydrogen 
break points, p curve, helium 
break points, p curve, nitrogen 
p_, absolute viscosity of gas, lb/(sec)(ft) 

S 
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PROGRAM LISTING 


- HEWMAN - 0003 

HEATED WIRE MESH ANALYSIS 0004 

BFTC HEWMAN HEWMOOOO 

HEWM000I 

ANALYTICAL INVESTIGATION OF ELECTRICALLY HEATED WIRE MESH FOR FL0WHEWM0002 
ING GAS HEATERS HEWM0003 

HEWM0004 

DIMENSION A (20.21) , ANSI 20) ,D135) , ASF (35) , DEI 35) ,A1 ( 2 ) * A2 ( 2 ) ,A3I2) ,HEWM0005 
1T2(2) ,TS(2) ,0(2) , RHO ( 2 ) , TB ( 2 ) .P(35) , TTS ( 35 ) , TT2 ( 35 ) ,QQQ(35) , DRES ( 3HEWM0006 
25) ,QIN(35) ,XLA(35) ,E(35) ,XIT(35) , IT (35) .TITLE (6) .TEST (3) .CPBPH13) .HEWM0007 
3CPAH13.5) .XKBPHI4) , XK AH ( 4, 5 ) , XMBPH ( 2 ) , XMAH ( 2 , 5 ) , C PBPHE ( 1) .CPAHE ( 1 .HEWM0008 
45) ,XKBPHE(5),XKAHE(5,5),XMBPHE(5) , XMAHE < 5, 5 ) , CPBPN < 3 > , C P AN ( 3 , 5 ) , XK.HEWM0009 


5BPN(3),XKAN(3,5) , XMBPN ( 2 ) , XMAN ( 2 , 5 ) ,CPA(5) , XKA ( 5 ) , XMA ( 5 ) , OG ( 2 ) HEWM0010 
DOUBLE PRECISION A.ANS HEWM0011 

DATA T 1TLE/6HHYDR0G,3HEN .6HHELIUM.3H , 6HN I TROG , 3HEN / HEWM0012 

DATA TEST/2HH2.2HHE.2HN2/ HEWM0013 

HEWM0014 

CURVE FIT COEFFICIENTS HEWM0015 

HEWMOOlfe 


DAT A ( CPBPH ( I ) , ( C P AH ( I ,J ) ,J = 1,5) , 1 = 1 ,3 1/1600., . 60895594 E- 13, -.17668HEWM0017 
1624E-09..19863776E-06.-.58512098E-04, . 34488570E+01 , 3400 . , .667 1 1 35 1HEWM00 18 
2E-14.-.96354453E-10, .47 1 86254E-06 , 68641 899E-03 , . 37794389E+01 , 999HEWM001 9 

39. . .33788 12 7E-14,-. 5440 1793E-10, .29369051 E-06,-. 4062 7483E-03, . 3680HEWM0020 

43507E+01 / HEWM002 1 

DATA ( XKBPH ( I ) , ( XKAH ( I , J ) , J=1 , 5 ) , I =1 , 4 ) / 3240 . , -.44762447E-1 8 , .45602HEWM0022 
1272E-14.-.14757960E-10, .4975 8692E-07 , . 67895535E-05 , 4200 . ,- . 2244840HEWMI0023 
24E-17, .19534089E-1 3, -.140343 18 E-10,-. 2025241 2E-06, .506367 88E-03 , 50HEWM0024 

300. . . 1 342 1 374E-1 7, -. 79444732 E-14,-. 1 5509958E- 1 0 , . 2 1 1 92 557E-06 . 28HEWM002 5 

4865473E-03.9999, , . 3 1 1 5352 8E-18 , . 142547 1 1 E-14 ,-. 79974343 E-l 1 ,-. 5965HEWM0026 
59474E-07, . 3 5433426E-03/ HEWM0027 

DATA (XMBPH (I) , (XMAH II, J) ,J=1,5) , I =1 , 2 ) / 3240 . , -. 412 57712E- 19 , .46966HEWM0028 
1670E-15.-.21505826E-11 , . 9 1476697E-08 , . 16277624E-05, 9999. , . 16330742HEWM0029 
2E-19,-.23379722E-15,.12630465E-11, . 13345724E-08 , . 868941 17E-05/ HEWM0030 

DATA(CPBPHE( I ) , (CPAHEI I ,J ) ,J=1,5) , 1=1, 11/9999, ,0.0,0.0,0.0,0.0,1.2HEWM0031 
148/ HEWM0032 

DATA ( XKBPHE ( I ) , (XKAHEI I , J ) » J=1 ,5 ) , I = 1 , 5 ) / 1080 . , -. 77890956E-1 7 , . 277HEWM0033 
180897E-13,-.42174067E-10, . 56632 1 75E-07 , . 35 1 87003E-05 , 1 980 . ,-. 10875HEWM0034 
2264E-18, . 14469 52 8 E-14,-. 7591 7067E-1 1 , . 35675647E-07 , . 85407794E-05 , 3HEWM0035 

3060.. . 22 11 5952 E-l 8,-. 2009 1672E-14,. 543 16253E-11,. 143839 10E-07, .213HEWM0036 

499514E-04.4320. , .9937 5672E-1 8, -.1440978 IE-1 3, .77 161 781 E-10,-. 16520HEWM0037 
5692 E-06, .18686650E-03.9999., . 10734377 E-l 8 ,-. 1 8710445E-14 , . 1 1455015HEWM0038 
6E-10.-.14133485E-07, . 5828377 5E-04/ HEWM0039 

DATA ( XMBPHE ( I ) , (XMAHE (I, J) , J = l,5) , I = 1 , 5 ) / 1 080 . , - .43709039E-1 7 , . 1 54HEWM0040 
1 62 508E-13.-.23183404E-10, .30621 850E-07, . 1 8674307 E-05 , 1980. ,- .49656HEWM0041 
2 867 E-l 8,. 3 4463889 E-l 4, -.10082 049 E-10,. 2 5066486 E -07,. 245 087 96 E-0 5, 3HEWM0042 

3060. . .96 14006 7 E-l 9,-. 8079 2 8 86E-15, . 1 7373264E- 1 1 , . 9949 1 83 1 E-0 8 , . 996HEWM0043 
406022E-05, 4320., .691 56058E-1 9, -. 97000761 E-l 5, .465 64625E-11,-. 37 171 HEWM0044 
582 5E-09, .212 18701 E-04, 9999. , . 4548 0686 E- 19 ,-. 9089472 8 E-l 5 , . 65060669HEWM0045 



ooo non ooo 


6E-ll,-.12201804E-07,.41129947E-04/ HEWM0046 

DAT A ( CPBPN { I ) , (CPAN( I , J ) , J = 1 , 5 ) , I = 1 , 3 ) / 1 200. , - . 32947 199E-13 , . 1 1 359HEWM0047 
1297E-09,-.11757656E-06, .51070 120E-04, .24025236,2600. ,- .38860908E-1HEWM0048 

24. . 2 7320 397E- 10,— *76523934E-07 , * 12770844E-03, • 1781 1544,9999* , .3465HEWM0049 
39852E-15,-.52774543E-11, . 275277796-07, -.49201293E-04, .31484032/ HEWM0050 

DATA ( XKBPN ( I ) , ( XK AN ( I , J ) , J = 1 , 5 ) , I = 1 , 3 ) / 1 900. , - . 42 1 7336 8E- 1 8 , . 22 349HEWM005 1 
1799E— 14,-.44977566E-1 1, .9042 1007E-08 , . 3570 8245E-06, 3800 . , .48366958HEWM0052 
2 E-l 9,-. 55485031 E-l 5, .203062 18E-11 , .16420 137E-08, . 38586614E-05 , 9999HEWM0053 

3., -.58176164E-19, . 1 0634372E-14, -.73856933 E-l 1 , . 26436289E-07 ,-.2095HEWM0054 

49671E-04/ HEWM005 5 

DATA ( XMBPN ( I ) , ( XMAN ( I , J ) , J=1 , 5 ) ,1=1,2) 71800* . 962 13327E-1 8 , . 52937HEWM0056 
1485 E— 14, -.12553952 E- 10,. 252081 25 E-07 ,.11203481 E-05 , 9999. , - . 1 637768HEWM0057 


23E-1 9, . 29989957E-1 5, -.23507248 E-l 1, .1 582500 6 E-07, .41787548E-05/ 
READ INPUT 

1 READ (5,107) CODE , NN 

READ (5,100) (D(I),E(I),ASF(I),DE(I),P{r),XLA(I),XIT(I),I=l,NN) 
READ (5,100) PCE,AR,XL,FR,T11 

PRINT TITLE 

1=0 

DO 8 1 = 1, 6, 2 
L = L+ 1 

I E ( CODE .NE . TEST ( L ) ) GO TO 8 
11=1+1 

WRITE (6,101) (TITLE ( J ) , J=I , II ) 

III=II/2 
8 CONTINUE 

STARTING VALUES FOR NEWTON-R APHSON 

TSG=1 500. 

T 1G=T 1 1 
1 1 1 1 = 1 

97 RH0G=5.2E-6*(TSG*1.0E-3)**1.185 

QG( 1 )=( E( 1 )**2*.948E-3)/(RH0G*XLA( 1 )*XL) 

GO TO <90, 91, 92), III 

90 CPG=3 .47 
GO TO 93 

91 CPG=1 .248 
GO TO 93 

92 CPG=. 248 

93 T2G=(QG( 1 ) +C PG*T 1G*FR ) / ( FR*CPG ) 

TBG=(TlG+T2G)/2.0 

CTS=TSG 
GO TO 94 

95 XG1=(FR*144.0*DE(1 ) ) / ( P ( 1 )*AR*XM ) 

XG2=( AR*ASF( 1 )*.364)/ ( 144.0*DE( 1 ) ) 


HEWM0058 
HE WM0059 
HEWM0060 
HEWM0061 
HEWM0062 
HEWM0063 
HEWM0064 
HEWM0065 
HEWM0066 
HE WM0067 
HEWM0068 
HEWM0069 
HEWM0070 
HEWM007 1 
HEWM0072 
HEWM0073 
HEWM0074 
HEWM007 5 
HEWM0076 
HEWM0077 
HEWM0078 
HEWM0079 
HEWM0080 
HE WM008 1 
HEWM0082 
HEWM0083 
HEWM0Q84 
HE WM0085 
HEWM0086 
HEWM0087 
HEWM0088 
HE WM0089 
HEWM0090 
HEWM0091 
HEWM0092 
HEWM0093 
HEWM0094 
HEWM0095 


CO 



r> r> o o o o 


00 

CO 



QG(2)=XG2*XK*(TSG-TBG)*XG1**.542 

HEWM0096 


I F ( OG ( 2 ) • GT * QG ( 1 ) ) GO TO 96 

HEWM0097 


TSG=TSG+ 100 • 

HEWM0098 


IF(TSG.LT*5800. ) GO TO 97 

HEWM0099 


T SG= 1 500 • 

HEWM0100 


T2G=600. 

HEWM0101 


QG(2)=1.0 

HEWM0102 

96 

TS I N=TSG 

HEWM0103 


T 2 I N=T2G 

HEWM0104 


II I 1=0 

HEWM0105 


DO 98 1=1, NN 

HEWM0106 


Q I N ( I ) =QG ( 2 ) 

HEWM0107 

98 

CONTINUE 

HEWM0108 


NNTEST=XIT (NN ) 

HEWM0109 


RES=0 • 0 

HEWM01 10 


T21=T2IN 

HE WMO 111 


TS1=TS IN 

HEWM0112 


DO 3 1=1, NN 

HEWM0113 


QOQ ( I ) =Q I N ( I ) 

HEWM0114 

3 

CONTINUE 

HEWM0115 


KKC = 1 

HEWM01 16 


J = 1 

HEWM01 17 
HEWM01 18 


PRINT INPUT FOR EACH STAGE 

HEWM0119 

HEWM0120 

5 

WRITE (6,102) KKC 

HEWM0121 


WRITE (6,106) 

HEWM0122 


WRITE (6,103) 

HEWM0123 


WRITE (6,104) D( J ) , E ( J ) , ASF ( J ) , DE ( J ) , XL , P ( J ) , AR 

HEWM0124 


WRITE (6,109) 

HEWM0125 


WRITE (6,104) FR,PCE,T11,XLA(J) ,T21,TS1 ,QQQ( J ) 

HEWM0126 


WRITE (6,108) 

HEWM0127 


T 1 =T 1 1 

HEWM0128 


T2 ( 1 ) = T21 

HEWM0129 


TS ( 1 ) = TS1 

HEWM01 30 


0(1)= OQQ(J) 

HEWM0131 


KTR = 1 

HEWM0132 

19 

K=1 

HEWM0133 

HEWM0134 


PROPERTIES OF ELEMENTS COMPUTED 

HEWM01 35 
HEWM0136 


GO TO (40, 41, 42), III 

HEWM0137 

40 

DO 43 11=1,3 

HEWM0138 


IF(T1.GT.CPBPH( 11 ) )G0 TO 43 

HEWM0139 


DO 44 Jl=l , 5 

HEWM0140 


CPA(J1)=CPAH( I 1 , J 1 ) 

HEWM0141 

44 

CONTINUE 

HEWM0142 


GO TO 49 

HEWM0143 

43 

CONTINUE 

HEWM0144 


WRITE (6,114) T 1 

HEWM0145 



GO TO 1 

41 DO 45 J 1 = 1 , 5 

1F(T1.GT.CPBPHE(1>) GO TO 48 
CPA(J1)=CPAHE(1,J1) 

45 CONTINUE 
GO TO 49 

48 WRITE (6fll5)Tl 


HEWM0146 

HEWM0147 

HEWM0148 

HEWM0149 

HEWM0150 

HEWM0151 

HEWM0152 


GO TO 1 

42 DO 46 11=1,3 

I F ( T1 .GT.CPBPN ( II ) )GO TO 46 

DO 47 J 1 =1 , 5 

CPA ( J 1 ) =C PAN ( 1 1 , J 1 ) 

47 CONTINUE 
GO TO 49 
46 CONTINUE 

WRITE (6,116) T 1 

49 CP=0.0 

DO 50 11=1,5 
CP=CP*T1+CPA( II ) 

50 CONTINUE 
CTS = T S { K ) 

94 GO TO (51, 52, 53), III 

51 DO 55 11=1,4 
IF(CTS.GT.XKBPH( II ) ) GO TO 55 
DO 56 J 1 =1 , 5 

XKA ( J 1 ) =XKAH ( 1 1 , J 1 ) 

56 CONTINUE 
GO TO 59 

55 CONTINUE 

WRITE(6,117) CTS 
GO TO 1 

52 DO 57 11=1,5 

IF(CTS.GT.XKBPHE( il) ) GO TO 57 

DO 58 J 1=1 , 5 

XKA ( J I )=XKAHE( 1 1 , J 1 ) 

58 CONTINUE 
GO TO 59 

57 CONTINUE 

WRITE (6,118) CTS 
GO TO 1 

53 DO 60 11=1,3 
IF(CTS«GT.XKBPN( II) ) GO TO 60 
DO 61 J 1=1 , 5 
XKA(J1)=XKAN( II ,J1 ) 

61 CONTINUE 
GO TO 59 
60 CONTINUE 

WRITE(6,119) CTS 
GO TO 1 

59 XK=0 • 0 


HEWM0153 
HEWM0154 
HEWM0155 
HEWMOl 56 
HEWM0157 
HEWMOl 58 
HEWM0159 
HEWM0160 
HEWM0161 
HEWMOl 62 
HEWM0163 
HEWM0164 
HEWM0165 
HEWM0166 
HEWM0167 
HEWM0168 
HEWM0169 
HEWM0170 
HEWM0171 
HEWMOl 72 
HEWM0173 
HEWM0174 
HEWM0175 
HEWMOl 76 
HEWM0177 
HEWMOl 78 
HEWM0179 
HEWM0180 
HEWM0181 
HEWM0182 
HEWM0183 
HEWMOl 84 
HEWM0185 
HEWM0186 
HEWM0187 
HEWM0188 
HEWM0189 
HEWM0190 
HEWM0191 
HEWM0192 
HEWM0193 
HEWMOl 94 
HEWM0195 



O o o 


CO 


DK=0 • 0 

DO 62 I 1=1,5 
XK=XK*CTS+XKA( I I ) 

IHII.EO.5) GO TO 62 
DK=DK*CTS+XKA{ II ) *FLOAT ( 5- 1 I ) 

62 CONTINUE 

GO TO (63,64,65) y I 1 1 

63 DO 66 11=1,2 
IF(CTS.GT.XMBPH(I1)) GO TO 66 
DO 67 J 1 = 1 , 5 

XMA ( Jl )=XMAH( 1 1 ,J1 ) 

67 CONTINUE 
GO TO 69 

66 CONTINUE 

WR I TE ( 6 , 1 20 ) CTS 
GO TO 1 

64 DO 68 11=1,5 

IF(CTS.GT.XMBPHE( ID ) GO TO 68 

DO 70 J 1=1 , 5 

XMA< Jl)=XMAHE< II, Jl) 

70 CONTINUE 
GO TO 69 

68 CONTINUE 

WR I TE ( 6 , 121 ) CTS 
GO TO 1 

65 DO 71 11=1,2 
IHCTS.GT.XMBPNI II ) ) GO TO 71 
DO 72 J 1 =1 , 5 

XMA(J1 )=XMAN( II, Jl ) 

72 CONTINUE 
GO TO 69 

71 CONTINUE 

WRI TE ( 6 , 1 22 ) CTS 
GO TO 1 

69 XM=0 mO 
DM=0 *0 

DO 73 11=1,5 
XM=XM*CTS+XMA< I I ) 

I F ( I I * EQ * 5 ) GO TO 73 
DM=DM*CTS+XMA( I I ) AFLOAT ( 5-1 1 ) 

73 CONTINUE 

I F ( I I I I • EQ • 1 ) GO TO 95 
IF(TS(K) .GT.O.O) GO TO 74 
WRITE(6,123) TS { K ) 

GO TO 1 

SET UP MATRIX OF COEFFICIENTS 

74 RHO(K) = 5*2E-6*(TS(K)/1000.0)**1.185 

A1(K) = ( ( E ( J )**2 > ! S «948E~3 ) / ( RHO (K)*XLA(J)#XL ) ) — Q ( K ) 


HEWM0196 

HEWM0197 

HEWM0198 

HEWM01 99 

HEWM0200 

HFWM0201 

HEWM0202 

HEWM02 03 

HEWM0204 

HEWM0205 

HEWM0206 

HEWM0207 

HEWM0208 

HEWM0209 

HEWM0210 

HEWM02 1 1 

HEWM0212 

HEWM02 1 3 

HEWM02 14 

HEWM02 1 5 

HEWM0216 

HEWM02 17 

HEWM0218 

HEWM02 19 

HEWM0220 

HEWM0221 

HEWM0222 

HEWM022 3 

HEWM0224 

HEWM022 5 

HEWM02 26 

HEWM022 7 

HEWM0228 

HEWM0229 

HEWM0230 

HEWM0231 

HEWM0232 

HEWM0233 

HEWM0234 

HEWM02 35 

HEWM0236 

HEWM0237 

HEWM0238 

HEWM02 39 

HEWM0240 

HEWM0241 

HEWM0242 

HEWM0243 

HEWM0244 

HEWM0245 



30 


A 2 ( K ) = FR*CP*(T2(K)-T1)-Q(K) 

TB ( K ) = <T1+T2(K) )/2.0 
XC0N1=(FR*144.*DE( J ) )/(P( J)*AR*XM) 

XC0N2=( AR*ASF( J )*.364)/< 144.*DE( J ) ) 

A3(K) = (XC0N2#XK 5 MTS(K)-TB(K)) 5 4 £ XC0N1#*.542)-Q(K) 
DA=XC0N2*(XK+DK*(TS(K)-TB(K) ) ) 

DB=( .542/ (XC0N1**.458) )*( < -XC0N 1*DM ) / XM ) 

DO 30 11=1,3 
DO 30 M=l,4 
A ( I I , M) =0 .ODO 
CONTINUE 


A ( 1 , 1 ) 

= 

- 1.0 

A ( 2 , 1 ) 

= 

- 1.0 

A ( 3, 1 > 

= 

- 1.0 

A(l, 2) 

= 

<-1.185*E(. 

7(5. 2E- 

-6#XLA(J )*XL) 

A { 2 ,2 > 

= 

0.0 

A ( 3 , 2 ) 

= 

XC0N2*XK*( 

A(l, 3) 

= 

0.0 

A ( 2 , 3 ) 

= 

FR*CP 

A ( 3 , 3 ) 

= 

- • 5*XC0N2* 

Ad, 4) 

= 

-A1 (K) 

A ( 2 ,4 ) 

= 

— A2 ( K ) 

A ( 3,4) 

= 

—A3 ( K ) 


>*DB+XC0N1**.542*DA 


14 

16 

13 


18 


CALL MGAUSD ( A 1 3 9 ANS ) 

AH=ANS ( 1 ) 

AK = ANS ( 2 ) 

AL = ANS ( 3 ) 

KP=K+1 

0 ( KP ) =0 ( K ) +AH 
TS ( KP ) =TS ( K )+AK 
T2 ( KP ) =T2 ( K ) +AL 
TEST0=ABS((Q(l)-0(2))/Q(2)) 

TESTTS = ABS( ( T S ( 1 ) -TS ( 2 ) )/TS(2) ) 
TESTT2=ABS ( ( T2 { 1 ) -T2 ( 2 ) ) / T2 ( 2 ) ) 
KTR=KTR+1 

IF(KTR.GT.SO) GO TO 28 
IF(TESTO-PCE) 14,13,13 
IF(TESTTS-PCE) 16,13,13 
I F ( TESTT2-PC E ) 18,13,13 
0 ( 1 ) =0 ( 2 ) 

TS(1)=TS(2) 

T2 ( 1 ) =T2 ( 2 ) 

GO TO 19 
TTS ( J ) =TS ( 2 ) 

TT2 < J ) =T2 ( 2 ) 

OOQ ( J ) =0 ( 2 ) 

RH0(2)=5.2E-6*(TS(2)/1000.0)**1.185 
DRESl J >=1.0/ (XLA( J)*RH0(2 )*XL) 
RES=RES+DRES(J ) 


HEWM0246 
HEWM0247 
HEWM0248 
HEWM0249 
HEWM02 50 
HEWM0251 
HEWM02 52 
HEWM0253 
HEWM0254 
HEWM0255 
HEWM02 56 
HEWM02 57 
HEWM0258 
HEWM0259 
►185 ) HEWM0260 
HEWM0261 
HEWM0262 
HEWM0263 
HEWM0264 
HEWM0265 
HEWM0266 
HEWM0267 
HEWM0268 
HEWM0269 
HEWM0270 
HEWM0271 
HEWM0272 
HEWM0273 
HEWM0274 
HEWM0275 
HEWM-0276 
HEWM0277 
HEWM0278 
HEWM0279 
HEWM0280 
HEWM0281 
HEWM0282 
HEWM0283 
HEWM0284 
HEWM0285 
HEWM0286 
HEWM0287 
HEWM0288 
HEWM0289 
HEWM0290 
HEWM0291 
HEWM0292 
HEWM0293 
HEWM0294 
HEWM02 95 



O o o 


CO 

o 


c 

C 

C 


HEWM0296 
HEWM0297 
HEWM0298 
HEWM0299 
HEWM0300 
HEWM030 1 
HEWM0302 
HEWM0303 
HEWM0304 
HEWM0305 
HEWM0306 
HEWM0307 
HEWM0308 
HEWM0309 
HEWM0310 
HEWM0311 
HEWM0312 
HEWM03 13 
HEWM0314 
HEWM0315 
HEWM0316 
HEWM03.17 
HEWM03 18 
HEWM0319 
HEWM0320 
HEWMQ32 1 
HEWM0322 
HEWM032 3 

100 FORMAT ( 7E8 * 5 ) HEWM0324 

101 FORMAT (84H1ANALYTICAL INVESTIGATION OF ELECTRICALLY HEATED WIRE MEHEWM0325 

1SH FOR FLOWING GAS HEATERS , , A6 , A3 , 10HPR0P ERT I ES ) HEWM0326 

102 F0RMAT(6H0STAGE,I3) HEWM0327 

103 F0RMAT(1H0,7X,1HD,15X,1HE,13X,4H ASF , 1 3X , 2HDE t 1 AX , 2HXL , 14X , 2H P,14HEWM0328 

IX, 2HAR ) HEWM032 9 

104 FORMAT { 8E 1 6 • 8 ) HEWM0330 

105 F0RMAT(1H0,7X,2HT2,14X,2HTS,14X,2H Q,14X,2HRT) HEWM0331 

106 FORMAT ( 6H0 INPUT ) HEWM0332 

107 FORMAT ( A2 * 1 2 ) HEWM0333 

108 FORMAT (7H00UTPUT) HEWM0334 

109 FORMAT ( 1H0,6X,2HFR, 1 3X f 3HPCE f 1 3X , 3HT 1 1 , 1 3X , 3HXL A f 1 2 X , 8HT2 GUESS , 7XHEWM0335 


RT=1. O/RES 

PRINT OUTPUT FOR EACH STAGE 
WRITE (6*105) 

WRITE (6,104) T2(2),TS<2),G(2),RT 
I F ( TTS ( J ) .LE.5800.0) GO TO 24 
26 WR I TE ( 6 * 1 1 1 ) 

GO TO 1 

24 I T { J ) =X I T ( J ) 

I F ( KKC .NE * I T ( J ) ) GO TO 7 
WRITE<6,112) I T < J ) 

WRITE(6,126) 

RES=0.0 

I F { I T ( J ) . EQ . NNTEST ) GO TO 1 
7 KKC = KKC + 1 
TS1 = TT S ( J ) 

T21=TT2 ( J )+(TT2 ( J )-Tl ) 

T11=TT2 ( J ) 

I F ( KKC • GT • I T ( J ) ) J = J + 1 
GO TO 5 

28 WRITE (6,113) 

WRITE (6, 124) TS( 1) ,TS(2) ,T2( 1) ,T2(2) ,0(1 ) ,0(2) 
WRITE(6,125) TESTTS,TESTT2,TEST0 
GO TO 1 

FORMAT STATEMENTS 


1,8HTS GUESS, 9X,7H0 GUESS) 

111 FORMAT ( 55H0TEMPERATURE EXCEEDS 5800 DEGREES RANKINE 

112 FORMAT ( 6H0X I T =,I3) 

113 FOR MAT(55H0 ITERATIONS EXCEED 50 

114 FORMAT ( 13H0ARGUMENT T1=,E16.8,4 
IFOR HYDROGEN) 

115 FORMAT ( 13H0ARGUMENT T1=,E16.8,A 
IFOR HELIUM ) 

116 FORMAT (13H0ARGUMENT T1=,E16.8,4 
IFOR NITROGEN) 


EXCEEDS 

LIMIT 

OF 

T 

VS. CP 

CURVE 

EXCEEDS 

LIMIT 

OF 

T 

VS. CP 

CURVE 

EXCEEDS 

LIMIT 

OF 

T 

VS. CP 

CURVE 


HE.WM0336 
) HEWM0337 
HEWM0338 
) HEWM0339 
HEWM0340 
HEWM0341 
HEWM0342 
HEWM0343 
HEWM0344 
HEWM0345 





117 

FORMAT (13H0ARGUMENT 
1 OR HYDROGEN ) 

TS=,E16.8,45H 

EXCEEDS 

LIMIT 

OF 

T 

VS. 

K 

CURVE 

FHEWM0346 

HEWM0347 

118 

FORMAT (13H0ARGUMENT 
1 OR HELIUM ) 

TS=,E16.8,45H 

EXCEEDS 

LIMIT 

OF 

T 

VS. 

K 

CURVE 

FHEWM0348 

HEWM0349 

119 

FORMAT (13H0ARGUMENT 
1 OR NITROGEN ) 

TS=,E16.8,45H 

EXCEEDS 

LIMIT 

OF 

T 

VS. 

K 

CURVE 

FHEWM0350 
HEWM035 1 

120 

FORMAT (13H0ARGUMENT 
IFOR HYDROGEN) 

TS=,E16.8,45H 

EXCEEDS 

LIMIT 

OF 

T 

VS. 

MU 

CURVE 

HEWM0352 

HEWM0353 

121 

FORMAT (13H0ARGUMENT 
IFOR HELIUM ) 

TS=,E16.8,45H 

EXCEEDS 

LIMIT 

OF 

T 

VS. 

MU 

CURVE 

HEWM0354 

HEWM0355 

122 

FORMAT ( 13H0ARGUMENT 
IFOR NITROGEN) 

TS=,E16.8,45H 

EXCEEDS 

LIMIT 

OF 

T 

VS. 

MU 

CURVE 

HEWM0356 

HEWM0357 

123 

FORMAT (36H0PR0GRAM CANNOT PROCEED 


TS =,l 

E16. 

8) 




HEWM0358 


124 FORMAT { 8H0TS ( 1 ) =»E14.6»8H TS(2) =,E14.6,8H T2(l) = , 

1 =,E14.6,7H 0(1) =,E14.6,7H 0(2) =,E14.6) 

125 FORMAT (9H0TESTTS =,E14.6,9H TESTT2 =,E14.6,8H TESTQ =,E14.6) 

126 FORMAT ( 1 HL ) 

END 

$ I BFTC MGAUS 

SUBROUTINE MGAUSD ( A , N , ANS ) 

C THIS SUBROUTINE SOLVES FROM 2-6 

C SIMULTANEOUS LINEAR EQUATIONS 

DIMENSION A(20,21) t ANS(20) 

DOUBLE PRECISION A , ANS 

DO 1 1=1, N 
1 ANS ( I ) = 0 * ODO 
DO 10 I=1,N 
DO 9 J = I , N 

A( I, J + 1)=A( I,J + 1)/A(I,I ) 

IF(I-N) 9,20,9 
9 CONTINUE 
K= 1 + 1 

DO 8 II=K,N 
DO 8 JJ=I,N 

A( II ,JJ+1)=-A< II , I )*A< I ,JJ+1)+A( II,JJ + 1 ) 

CONTINUE 


8 

10 

20 


30 

31 


12 


11 


ANS ( N ) =A ( I ,J+1 ) 

IF(N-1 ) 31,30,31 
RETURN 
J=N-1 
I I=J 

DO 11 1=1,11 
K=J + 1 

DO 12 M= 1 , I 

ANS( J)=ANS(K)*A( J,K)+ANS( J > 
K = K+ 1 

ANS(J)=A(J,K) -ANS ( J ) 

J = J-1 

RETURN 

END 


E14.6,8H T2 ( 2 ) HEWM03 59 
HEWM0360 
HEWM036 1 
HEWM0362 
HEWM03-63 
MGAUOOOO 
MGAUOOO 1 
MGAU0002 
MGAU0003 
MGAU0004 
MGAU0005 
MGAU0006 
MGAU0007 
MGAU0008 
MGAU0009 
MGAU0010 
MGAUOO 1 1 
MGAUOO 1 2 
MGAUOO 1 3 
MGAUOO 14 
MGAUOO 1 5 
MGAUOO 16 
MGAU0017 
MGAUOO 1 8 
MGAUOO 19 
MGAU0020 
MGAU002 1 
MGAU0022 
MGAU002 3 
MGAU0024 
MGAU002 5 
MGAU0026 
MGAU0027 
MGAU0028 
MGAU002 9 
MGAU0030 
MGAU003 1 



CO 

00 


ANALYTICAL INVESTIGATION OF ELECTRICALLY HEATED WIRE PESH FCR FLOWING GAS HEATERS, NITROGEN PROPERTIES 

STAGE 1 

INPUT 

D _E ASF CE XL P AR 

C .3C CCOOCCE-C I C.lrfS/CCCCE C2 C.552CGCCGE 01 0 • 59400000 E-02 0.10625000E Cl C«7C<CCCCCE CO C.4250000CE 01 

C.CSCCCCCCE CO C . 1 CCCCCOOE- Cl C.5CCCC0CQE 03 0.27400000E 03 0 .77715221E C2 C.24CCCCCCE 04 0«4 < i352i77E 02 

OUTPUT 

T l 1 C C RT 

0.7787901EE C3 C.33T7C069E C4 C.47809145E 02 0 . 640 3 1 472 E-0 2 

STAGE 2 
INPUT 

C E ASF CE XL P AR 

C .3CC00U00E-C 1 C.1/97CCCCE C2 G.552C0CCCE Cl 0 • 594000 00 E-02 0.10625000E Cl 0.7C4C0CCCE 00 C.425000C0E 01 

C.69CC00CCE CC C . 1CCCCCC0E-C1 C.77879018E 03 Q.27400000E 03 0.10§75804^ C* C .11 -J^Yllle C4 02 

OUTPLT 

0.1C441995E C4 0 .3499751 SE C4 C.45E28684E 02 0 . 326928 72 E-02 
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